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EXPLICIT VG GUIDANCE ALGORITHM FOR 
A SOLID POWERED CLOSED LOOP 
GUIDANCE MISSION 

$. Dasgupta and S, Ramakrishna 

Control and (inidanrr Division, Ythrum Sarabhat Sfnur Trivandrum, India 


Ab s tract . An explicit form of velocity to be gained guidance algorithm is 
proposed for ASLV mission of Indian Space Research Organisation, ASLV 
mission is to inject a 150 kg payload into a 400 Km near circular orbit 
us^ng a four stage solid propellant system* Since thrust termination 
cannot be commanded when the lower stages overperform during the flight* 
guidance algorithm commands the final stage to be reoriented* such that 
the final orbit is close to circular* When underperformance occurs, a 
velocity augmentation system is activated to pull up the performance to 
near nominal* It is shown through simulation that for 1 sigma over and 
under performance, the orbit eccentricity is less than 0*00055* 

Key words ; Guidance systems; closed loop systems; space vehicles; 
target tracking; required velocity, 
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Topo centric Inertial coordinates parallel to LPI, 
witn origin at Earth Centre, 

Requiredvelacity in guidance co-ordinates* 

Gravitational constant* 

Distance of vehicle CG from earth centre 

desired angular momentum during long coast 

desired apogee at end of long coast. 

local horizontal component of velocity required. 

direction cosine of normal to the current orbit plane* 

position vector in LPI 

velocitv vector in LPI 

angle between desired orbit plane and current orbit plane. 
Constants 

steering command rate 
velocity to be gained vector 

acceleration vector 

pitch attitude command 
yaw attitude command 


normalised velocity to be gained components 


(W > p) 
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velocity and flight path angle at any time during long coast, 

equivalent velocity increment due to fourth stage, when 
added impulsively, 

angle,from local horizontal, at which fourth stage is 
ignited 

velocity decrement due to gravity over the duration of 
fourth stage burning, acting along local vertical* 

vector sum of ( ^ ) and ( 5 } 180°) 

Resultant velocity vector sum of (Va T ,ft ?r j 

and tv,,*,) 


I SSI 
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INTRODUCTION 

Multi stage launch vehicles are steered in 
flight in order to reach desired altitude, 
velocity, flight path angle and flight 
azimuth at the end of thrust, so that the 
payload is inserted into the desired orbit. 
While payload optimisation and Mission 
profile are generated assuming a nominal 
vehicle performance, the inflight performance 
of the vehicle's propulsion system and 
disturbances due to aerodynamics can cause 
significant dispersions in the orbit 
attained. This calls for a closed loop 
guided system that can measure position, 
velocity and attitude of the vehicle 
at any time and determine from them, the 
steering profile to be followed in order to 
achieve the desired orbit. Closed loop 
guidance algorithm is invoked to derive 
the steering commands in pitch, yaw and 
roll motions, 

Indian Space Research Organisation plans 
to have its first closed lodp guided 
mission with the launch vehicle ASLV, 

This is designed to put 150 Kg, payload 
into a near earth circular orbit with 
specified tolerances in orbit eccentricity 
and inclination, aslv is an all solid 
powered launch vehicle and has no provision 
for thrust termination, 

iXiring the atmospheric phase, the vehicle 
executes a predetermined pitching sequence* 
Closed loop guidance is initiated at second 
stage,and continues until the end of third 
stage burnout. This is followed by a long 
coast during which the vehicle gain altitude. 
At the apogee of the coast phase the vehi¬ 
cle is spun, fourth stage separated and 
ignited to impart the necessary velocity 
increment to the payload* 

Several guidance algorithms are reported 
in literature (Tounsend; Battln 1982 Gage 
1968), Broadly, these can be categorised 
into explicit and implicit schemes. In 
explicit guidance, closed form analytical 
solution is obtained for a simplified 
TPBV problem. Several assumptions, such as 
those ef uniform gravity, constant thrust 
and propellant massflow rate are made 
in order to arrive at analytical solutions. 

In implicit schemes, required velocities 
are obtained from a nominal profile and are 
stored as polynomials of perturbation in 
position, time, velocity etc. about the 
nominal burnout state. Generally, explicit 
schemes require higher computation time, 
while implicit schemes require a larger 
memory, 

GUIDANCE ALGORITHM FOR ASLV 

In this paper, a guidance algorithm is 
proposed, that is an implicit form of c 
classical velocity to be gained (VG) 
technique* While steering is achieved 
using cross product steering co^and 
of VG guidance, the required velocity is 
computed explicitly, using closed form 
analytical solutions* The target is 
specified as a desired coast apogee, that 
must be reached with a specified angular 
momentum. Using classical keplerian laws, 
the velocity required to attain the specified 
terminal conditions froM any altitude, 
are analytically evaluated. This obviates 
the need for storing predetermined nominal 
profile, of required velocities, either as 


tables or as polynomials, as is generally 
adopted in classical VG techniques. Since 
closed form Keplerian solution is used to 
predict the path connecting the current 
position of vehicle to the desired terminal 
state, the TPBV problem formulation is not 
necessary. Consequently# the memory 
requirement as well as computational load 
are expected to be less than conventional 
implicit and explicit schemes,respectively. 

In ASLV the fourth stage is a solid motor 
that is spin stabilised prior to ignition. 
There is no provision for thrust termina¬ 
tion In any of the stages of ASLV, 
Consequently# when any or all of the 
stages overperform, i.e, impart higher 
velocity than is nominally expected, the 
resultant orbit will be quite eccentric. 

I n order to achieve a near circular orbit, 
prior to ignition the final stage is 
reoriented from the local horizontal 
attitude by an angle determined on-board. 
This angle Is evaluated as a function of 
overperformance by 1ower stages* 

In case the lower stages under perform, a 
velocity Augmented system {VAS) mounted 
axially on ASLV, is fired during the coast 
phase. The orientation of the vehicle and 
the duration of VAS firing are determined 
en—board, depending on the level of 
under performance. 

The capacity of VAS, to Impart velocity 
is limited to 25 m/s, study indicates 
that under dispersions, it may be necessary 
to augment velocity by as much as 53 m/s* 
Left uncorrected, this will again result 
in an eccentric orbit. Effect of removing 
VAS (that weighs nearly 42 Kg) from ASLV 
and depending solely on the reorientation 
algorithm to achieve low eccentricity, la 
studied and results are presented in this 
paper, 

C nor dj natea Systern 

Guidance coordinates (Xq Yg Zg) are 
illustrated in Fig.l* Inertial coordi¬ 
nates ( X£ Yj Zj) are Topocsntric and 
defined along the launch point vertical, 
down range and normal to down range 
respectively. 

Required Velocity 

Target conditions are the states attained 
by nominal ASLV vehicle at the apogee of 
the third stage coast and are obtained 
from a nominal, optimised open loop tra¬ 
jectory* Using Keplerian equations to 
describe nature of a body under Central 
force field, for any given altitude, the 
velocity required to attain the specified 
altitude and velocity at cost apogee are 
given in guidance coordinates as 

^ 1_ = >1 r i + *.( tO" - ±a)/ a *J~ h m 

n-'- r 

if = H/X. 

Anqle-to-oo 

The plane of orbit, if coast were to 
commence at the current state, is given 
by I = S x The angle between this plane 
and the desired orbit plane is defined as 
the Angle-to-gos, V e loclty component vl, 
normal to the orbit plant, is then^where 
K1 is a constant^ ^ - K±-$ 
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Since navigation Is performed In LPI 
coordinates (u > v* , va) are transformed 
to LPI and compared with current veloci¬ 
ties in LPI, to give the velocity to be 
gained vector in LPI coordinate system. 

Crossproduct Steering 

Battin (1982) has shown that if the 
vehicle is steered according to 

wc - xfv'&xva) 

where is the commanded rate, 

Vg id the velocity to be gained 

V& is the derivative of VCc 

the path followed by the vehicle is fuel 
optimum. 

An alternate scheme is to use , the 
acceleration vector, instead of VG so 
that „ — — N . 

*J C - K ( VG * a ) (Z) 


local vortical) altitude of aslv is not 
significantly affected by the 'fine trim* 
phase of VAS firing. 

The 3 cr under performance level of ASLV 
lower stages Is equivalent of 53 m/s at 
third stage burnout. However, VAS tank 
capacity is limited to a velocity 
augmentation of only 25 m/s m As may be 
expected, for underperformances demanding 
velocity corrections beyond 25 m/a f the 
eccentricity of the resultant orbit Is 
high* and beyond the specified 0,0025, 

In order to hold eccentricity to within 

specifications, under ± 3rievels of vehicle 
dispersions, VAS, which weighs 42 Kg, has 
been removed from the structural weight of 
third stage, in the simulation. This 
leaves a level of overperformance of 
nearly 70 m/s. at third stage bum out. 
Since, the + 34T dispersion is 53 m/s, 
sufficient margin is available in the 
vehicle to handle the expected dispersions, 
and velocity augmentation is not called for. 


Thus, by computing Vdr analytically and 
steering through croseproduct,pitch and 
yaw commands are generated with small 
computational load, small data storage 
and being an explicit algorithm, gives 
accurate injection conditions. 

The flow chart for steering command 
generation during thrust phase is shown 
In fig,2. It may be observed from Pig.2 
that required velocity is initially 
computed in guidance coordinates in terms 
of 1 in—plane 1 and 'out-of-plane 1 velocity 
components. These are transformed to 
inertial co-ordinate systems for calcu¬ 
lating velocity to be gained in inertial 
system. Steering rate is then computed 
through (2), Near the end of guided 
phase, steering angle Is computed 
according to 

B c - ,{31 

~ - sin l [ /a r*f] (4) 

Details of mathematical formulation are 
discussed by E^sgupta and Ramakrishna 
(19821 

Closed loop Guidance In Coast Phase 


In case the lower stages under perform, 
the apogee and anqular momentum at apoaee 
are enhanced by VAS* At the 

beginning of coast, whenever the predicted 
apogee is in error greater than a preset 
amount, VAS is fired, IXiring this period 
of VAS actuation, the vehicle is oriented 
according to equations (3) and (4), It 
is also ensured that the vehicle tracking 
error is within the stipulated control 
system dead-zone, prior to VAS actuation, 
when the vehicle's velocity vector is nearly 
locally horizontal, VAS actuation is based 
on predicted velocity errors rather than 
apogee. This phase is termed 'fine trim' 
of VAS and velocity error of less than 
0,1 m/a and 50m can be achieved at the 
end of VAS actuation* During this phase, 
the vehicle is oriented to the locally 
horizontal* Since the vehicles flight 
path angle is nearly 90°(measured from 


I n case the lower stages overperform, the 
apogee of the coast phase will be higher 
than that desired. If the final stage 
is fired at this apogee, the resultant 
orbit may be highly eccentric. Pig,3 
gives the velocity when tfce final 
stage is ignited at an angle o c to the 
local horizontal. In actual practice, 
the velocity addition during the stage 
vHIl be along an arc subtending an angle 
£^at earth centre. However, for 
simplicity an impulsive addition of velo¬ 
city is assumed during fourth stage. V4 
is the effective Impulsive velocity 
Increment, the value being determined 
from simulation of a nominal fourth 
stage trajectory. 

By forcing the following conditions: 

Magnitude of Vfe equal to the predicted 
circular 


velocity at end of final stage, and 

locally horizontal at burn out of 
final stage, where is the predicted 
velocity vector at fourth stage burnout. 


the reorientation angle oc and instant of 
final stage ignition is determined, such 
that the resultant orbit is nearly cir¬ 
cular, The final stage is reoriented and 
spun before ignition. In fig3, Sv^ is 
the decrement In velocity due to 
acceleration over the duration of final 
stage thrust, is the resultant of 

{V 5 j p ) and { 5 Tt)y^is estimated as 

*4 - o y £§T > ^ e re 

tf is the burning time of stage. 

At any instant, let R be the position of the 
vehicle. Circular velocity Vc at this 
altitude is V c - /M/si where y- is the 
gravitational constant. 


When V4 is added to V3 impulsively at an 
angle to local horizontal, such that 

~ ^ ^ )/z up*-/*] 

-■ (s) 

the resultant velocity is equal in magnitude 
to Vc, 

The flight path angle is given by 


ACS1905-C* 
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TA+C-c-V) 

where tan j- C v/ ?' 1 Us(£p.-t 

O.} 

In flight, the velocity addition is not 
impulsive but over the duration of fourth 
stage burning. It is estimated that the 
local vertical direction changes by 
where £*is nearly 2.3 deg over this 
duration. Thus in order to attain a 
flight path angle of 90° at fourth stage 
burn out, 

^ ” y ^ Sue, (g) 

During the coast phase, the vehicle pitch 
is continually commanded to the local 
horizontal* Equation (5) and (7) are solve! 
for**: and y respectively. Sz is evaluated 
at third stage, burnout as a function of 
residual velocity to be gained , Vjjbc 

(gj 

tfi t are precomputed through simulation, 
being 2,3 degress for ASLV, 

During coast phase, condition stipulated 
in equation (6) is checked at each guidance 
cycle, and when satisfied, the vehicle is 
reoriented through and coast phase is 
terminated, followed by spin up, separation 
and final stage ignition* Thus, it may be 
noted, that no prediction is made on the 
time at which fourth stage is to be 
ignited, 

SIMULATION RESULTS 

ASLV is designed to put 150 Kg payload in 
near 400 km, near circular orbit, eccen¬ 
tricity being less than . 002C The target 
for guidance algorithm is drawn from 
nominal fourth stage ignition states and 
ate 

Altitude 3400.644 KmS. 

Velocity =5198*1165 tv's, 

direction cosines 
of Normal to de¬ 
sired orbit plane^ =-0.0578, 0.99757, 
0.038696 

These are to be achieved at end of 3rd 
stage long coast. 

Table 2 gives the coast and condition 
reached through simulation, under various 
perturbations in thrust, weight and other 
parameters, RSS estimate of velocity 
error at 3rd stage burnout is 17,0 m/s. 

(1 sigma). For underperformances, VAS 
logic activates the axial thrusters, VAS 
trim is called when flight path angle is 
88 e . Prediced velo<fci£y error at the end 
of fine trim firing is 0,1 m/s. Altitude 
error is less than 50.0 m. 

The effect of out of plane disturbances 
due to thrust misalignment and winds is 
shown in Table 3. Under worst case of 
disturbances, the error in inclination is 
0.037* 


eccentricity of orbit for underperformances 
ata larger than 25 m/s, is beyond the 
specification of 0.0025. eg, for under¬ 
performance of 53 m/s, eccentricity is 
0*00749. Table 5 gives the orbital 
dispersions for ± 3<r dispersions, when 
VAS weight is removed from the third 
stage for the nominal vehicle- Vgbo=-7l, 2ny£ 
+3 o - dispersions in vehicle parameters 
correspond to a Vgbo of —122,9 m/s and 
-19.2 rr/s respectively. With reorientation 
angle oc indicated in Table 5, eccentricities 
are kept within the specification* 

Analytical expressions for sensitivity of 
Ap to errors in u,v are derived by 
Dasgupta (3), For ASLV, at third stage 
burn out, '3A J b/3*c =9 3,3 m/m/s* 

=74,8 m/m/s. Also 

a*vi*=_ # O033 m/s/m and =-.00082 

m/s/m. Thus 7 )a^a=-, 369226 

The navigation system of ASLV has an 
error bound of 300m at third stage 
burnout. Hence, the apogee reached at 
end of long coast is in error by a 
maximum of 110,76 m from the desired value* 

CONCLUSIONS 

This paper discusses the guidance stra¬ 
tegy that is followed for guiding an 
all-solid-propellant launch vehicle ASLV 
into a desired orbit. In addition to a 
steering during thrust phase, coast phase 
guidance with axial thruster firing and 
final stage reorientation Is proposed as 
the strategy to overcome the 
table effects of under and overperfor¬ 
mance of launch vehicle. Simulation 
results illustrate the performance of 
the algorithm. 
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When vehicle overperforms, the final stage 
is reoriented through an angle «c depen¬ 
ding on the level of overperformances* 

For a nominal final stage performance, 
the orbit reached for various levels of 
overperformances in lower stages Is given 
at 1 sigma dispersion in parameters, 
orbit eccentricity is 0,00054, When 

VAS capacity is limited to 25 rq/s. 
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TABLE 1 3 3lama PERTURBATION LEVELS 


no': Sta 9 e 

Parameter Disper¬ 
sion 

VG at 3rd 
stage burn¬ 
out (nv's) 

1 Strapons 

Thrust +1,0% 


-17,02 

2 Strapons 

Thrust -1*0% 


+13.86 

3,0ore 

Thrust 4-1,0% 


-19.34 

4, Core 

Thrust-1.0% 


+16.09 

5,Second 

Thrust +1.0% 


-19.66 

6, Second 

Thrust -1.0% 


+16.39 

7,Third 

Thrust +1.0% 


-20.46 

8.Third 

Thrust-1,0% 


+17.18 

9,Strapons 

Str.we.-20,0Kgs. 

- 3,40 

10,Strapons 

Str.wt.+20.0 

n 

- 0.06 

11. Core 

Str.wt.-20 

ii 

-5.32 

12 .Core 

str.wt.+20 

i< 

+1.97 

13.Second 

Str.wt.-15 

ii 

-7.56 

14 .Second 

str+ wt.+15 

n 

+4.19 

15,Third 

Str.wt-7,5 

ii 

-10.90 

16, Third 

Str*wt+7,5 

n 

+ 7.57 

17.Strapons 

Propt.wt.+30 

ii 

-5,85 

IS,Strapons 

Propt,wt,-30 

14 

+2.54 

35. Core 

Prop.wt. +30 * 


-4.98 


1 

2 

3 

4 

5 

13 

Second 

Str,wt. 

+1558,83 

+05,3435 

14 

Second 

Str.wt. 

-0774,53 

-03,1565 

15 

Third 

Str,wt. 

+2285,80 

+07,7135 

16 

Third 

Str.wt. 

-1389.77 

-05,7265 

17 

Strapons 

Prop.wt. 

+1190.22 

+04,1750 

18 

Strapons 

Prop.wt. 

-0460.39 

-01,8871 

19 

Core 

Froppwt, 

+1018.86 

+03.5411 

20 

Core 

Prop.wt, 

-0325,87 

-01.0220 

21 

Second 

Prop.wt, 

+1083.51 

+07.1330 

22 

Second 

Prop.wt, 

-1254.32 

-05,1426 

23 

THird 

Prop.wt, 

+2287,59 

+07,8473 

24 

Third 

Prop.wt, 

-1436,98 

-05.9033 

25 

First 

Drag.coeff+6150,09 

+21,1236 

26 

First 

Drag. “ 

-4880.30 

-20.0640 


TABLE 3 EFFECT OF YAW DISTURBANCE 

ON ORBITAL ERRORS 

Note; 1. Misalignments in thrust for first, 
second ind third stages are 0,3, 
0,28, 0,15 degrees respectively, 

2. Misalignment in yaw plane is 90° 


2D .Core 

Propt,wt.-3p n 

+1.65 

21 .Second 

Propt,wt,+20 " 

-10.10 

22 .Second 

Propt.wt.-20 " 

+ 6.81 

23.Third 

Propt.wt.+7,5 a = 

-11.10 

24.Third 

Propt.wt.-7.5 11 

+ 7.81 

25,First 

Drag.coeff -15*0% 

-30.06 

26.First 

Drag.coef f+15 11 

+26.53 


Yaw 

S* Dist. 

VT occurs 
No: . 


Incli¬ 

nation 

error 

(deg) 


Errors at injection 
(predicted) 


Apogee vel.at apogee 
(mtr) (m/a) 

(RT-AP) (VT-VAP) 


Re¬ 

marks 



1 . 

Strapons 

stage 

0,0003 

*71,196 

0,018 

RT, VT- 
target 
value AP, 

TABLE 2 INJECTION ERRORS FOR 3 SIGMA 

PERTURBATIONS 

2 , 

Core 

stage 

0,021 

829.319 

3,185 

VAP pre- 
dicted 

3. 

Second 

stage 

0,012 

-32.043 - 

0.310 

values. 

Injection errors 

SI, Stage Para- without cut off 

4. 

Third 

stage 

0,002 - 

-289.124 -1.305 


meter Apogee(mtrs) Vel.at . t 

5. 

All stage® 0. 037 

1089.629 

5.592 



1 

Strapons 

Thrust 

+0493,92 

+$2*0615 

2 

Strapons 

Thrust 

-2555.59 

-10,4765 

3 

Core 

Thrust 

+3974,50 

+13.6615 

4 

Core 

Thrust 

-2942.53 

-12.0613 

5 

Second 

Thrust 

+4130.97 

+13,7235 

6 

Second 

Thrust 

-2984.93 

-09.2365 

7 

Third 

Thrust 

+4286.83 

+14.3235 

8 

Third 

Thrust 

-3109,09 

-13.0865 

9 

Strapons 

Str.wt, 

+0703.13 

+02,3855 

10 Strapo rs 

Str.wt, 

+0016 ,11 

+00.0345 

11 

Core 

Str.wt. 

+1077.87 

+03.8035 

12Core 

Str.wt* 

-037,93 

-01.4565 
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TABLE 4 EFFECT OF OVERPERFORMANCE 

ON ORBITAL ERRORS 


si. 

Not 

VG 

(nv/s) 

Degrees 

Apogee 

Kms 

Perigee 

Kms 

Mean Dispersion 
Kins, Kms. 

Eccentricity 

1 

-lo.e 

6.930 

403.423 

410.090 

402.250 

+1.16 

0.00017 

2 

-16.8 

8.390 

406.032 

398.729 

402.300 

+3.60 

0.00054 

3 

-26.2 

10.166 

408.454 

394.727 

401.590 

+6.86 

0.00101 

4 

-39.3 

12.279 

410.928 

386.407 

398.677 

+12.26 

0.00181 

5 

-52.5 

14.135 

407.726 

378.383 

393.083 

+14.70 

0.00217 



ORBIT WITH NO 

TABLE 5 

VAS ON ASLV 

Vgbo m/s * 

deg 

Eccentricity 

-19.2 

9 .7 

0.0012 

-38.6 

12.3 

0.00023 

-53.12 

14.11 

0.00046 

-71.2 

15.9 

0.00133 

-89.21 

18.2 

0.00040 

-103.5 

19.0 

0.00146 

-122.9 

20.3 

0.00221 
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